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The infrared eye of the Wide-Field Camera 3 on the Hubble Space 
Telescope reveals multiple main sequences of very low-mass stars in 

NGC 2808. 1 

A. P. Milone 2 ' 3 , A. F. Marino 4 , S. Cassisi 5 , G. Piotto 6 - 7 , L. R. Bedin 7 , J. Anderson 8 , F. 
Allard 9 , A. Aparicio 2 ' 3 , A. Bellini 8 , R. Buonanno 5,10 , M. Monelli 2 ' 3 , A. Pietrinferni 5 

ABSTRACT 



^ . We use images taken with the infrared channel of the Wide Field Camera 

3 on the Hubble Space Telescope (HST) to study the multiple main sequences 
(MSs) of NGC 2808. Below the turn off, the red, the middle, and the blue MS, 
previously detected from visual-band photometry, are visible over an interval of 
about 3.5 F160W magnitudes. The three MSs merge together at the level of 
the MS bend. At fainter magnitudes, the MS again splits into two components 



O . 

containing ~65% and ~35% of stars, with the most-populated MS being the 
bluest one. Theoretical isochrones suggest that the latter is connected to the 
red MS discovered in the optical color-magnitude diagram (CMD), and hence 
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corresponds to the first stellar generation, having primordial helium and enhanced 
carbon and oxygen abundances. The less-populated MS in the faint part of the 
near-IR CMD is helium-rich and poor in carbon and oxygen, and it can be 
associated with the middle and the blue MS of the optical CMD. The finding 
that the photometric signature of abundance anticorrelation are also present in 
fully convective MS stars reinforces the inference that they have a primordial 
origin. 

Subject headings: globular clusters: individual (NGC 2808) - - Hertzsprung- 
Russell diagram 



1. Introduction 

In recent years, photometric studies have shown that the color-magnitude diagrams 
(CMDs) of globular clusters (GCs) can be very complex, with the presence of multiple main 
sequences (MSs, e. g. Anderson 1997, Bedin et al. 2004, Piotto et al. 2007, hereafter P07), 
multiple sub-giant branches (SGBs, e. g. Milone et al. 2008 and 2012a, Anderson et al. 2009, 
Piotto et al. 2012), and multiple or spread red-giant branches (RGBs, e. g. Yong et al. 2008, 
Marino et al. 2008, Lee et al. 2011). 

Photometric and spectroscopic investigations have revealed that the multiple sequences 
in the CMDs of many GCs are populated by stars with different helium and light-element 
abundances (e. g. Piotto et al. 2005, 2007, Yong et al. 2008, Marino et al. 2008, Sbordone 
et al. 2011, hereafter Sll). Is some cases, the presence of stellar populations with different 
composition (in particular, Helium) has been associated with the presence of multimodal 
or extended horizontal-branches (HBs) (e. g. D'Antona et al. 2005, P07, D'Antona & Caloi 
2008, Marino et al. 2011). 

Among clusters with multiple stellar populations, NGC 2808 is certainly one of the 
most intriguing objects. Its CMD shows a multimodal MS (D'Antona et al. 2005) composed 
of three distinct components (P07, Milone et al. 2012b, hereafter M12), a multimodal HB, 
which is greatly extended blueward (Sosin et al. 1997, Bedin et al. 2000), and a spread RGB 
(Lee et al. 2011). Furthermore, spectroscopic studies of RGB, HB, and bright MS stars 
have revealed significant star-to-star variations in the light-element abundances, with three 



1 Bascd on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope 
Science Institute, which is operated by AURA, Inc., under NASA contract NAS 5-26555. 
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distinct groups of stars populating an extended Na-0 anti-correlation (Carretta et al. 2006, 
Gratton et al. 2011, Bragaglia et al. 2010). 

Photometric studies, based on data collected with the Wide Field Channel of the Ad- 
vanced Camera for Survey (WFC/ACS) on board HST have made it possible to detect and 
characterize the multiple MSs of NGC 2808 from the MSTO down to about 4 magnitudes 
below the MSTO (P07, M12): the stars all have almost the same age and [Fe/H] but different 
Helium and light-element abundances. The red MS (rMS) corresponds to a first generation 
and has primordial helium and light-element abundances (Y~0.25), while the blue (bMS) 
and the middle (mMS) include later generations of stars, and are both enhanced in helium, 
sodium, and nitrogen (Y~0.38 and Y~0.32, respectively) and depleted in oxygen and carbon 
(D'Antona et al. 2005, P07). 

Usually, photometry of GC sequences extends over a limited spectral region, from the 
ultra-violet (A ~ 2000A) to the near- infrared (NIR, A ~ 8000A). As such, multiple sequences 
are rarely detected along the lower part of the MS, because observational limits make it hard 
to get high-accuracy photometry of very faint and red stars in optical and UV colors. 

In this Letter we use HST to extend the study to the near-infrared passbands and 
complement the work by Sll who analysed the signatures of the anticorrelations in the visual 
and ultra-violect portion of the spectrum. We analyze the CMD of NGC 2808 through the 
F110W and F160W (~i7) filters of the infra-red channel of the Wide Field Camera 

3 (WFC3/IR) and follow, for the first time, the multiple sequences of this cluster over a 
wide interval of stellar masses, from the turn off down to very low-mass (VLM) MS stars 
(M ~0.2M Q ). 

2. Observations and data reduction 

In this work we used archival HST images taken with the WFC3 /IR camera for program 
GO-11665 (P. I. Brown). This data set consists of 2x699s exposures through F110W and 
two through F160W of 799s and 899s for each of three different WFC3/NIR fields, taken 
in parallel while STIS was taking spectra of EHB stars. All fields are located at about six 
arcmin from the cluster center. 

In the following, fields 'A', 'B', and 'C will correspond, respectively, to the fields in 
the: south-west, south, and north-east of the cluster center. Field 'A' partially overlaps the 
ACS/WFC field analyzed by P07 and M12. Therefore, for a sub-sample of the measured 
stars, we have both visual and NIR photometry and can perform a direct comparison of the 
multi-band photometry (see Sect. [3]). 
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The images were reduced by using a software package that is based largely on the 
algorithms described by Anderson & King (2006) and will be presented in a separate paper. 
Star positions are corrected for geometric distortion by using the solution given by Anderson 
et al. (in preparation), and were calibrated as in Bedin et al. (2005) using the current on-line 
estimates for zero points and encircled energies^]. 

The analysis we present here requires high-precision photometry, so we selected a high- 
quality sample of stars that (1) have a good fit to the PSF, (2) are relatively isolated, (3) and 
have small astrometric and photometric errors (see Milone et al. 2009, Sect. 2.1 for details 
of this procedure). Finally we corrected our photometry for differential reddening using the 
method that is described in great detail in Milone et al. (2012c). Briefly, the method consists 
in defining a MS ridge line in the CMD. Then, for each star, we selected the 35 nearest well- 
measured neighbors and evaluated for each of them the color distance from the fiducial line 
along the reddening line. We applied to the target star a correction equal to the median 
color distance of these 35 stars. 



3. The NIR Color-Magnitude diagram of NGC 2808 

The IR m F110W versus mpnow ~ ^Fi60W CMDs of stars in the fields 'A', 'B' and 'C 
are plotted in the upper panels of Fig. [TJ while the lower-left and lower-right panels show 
the mFi6ow versus mFiiow — ^Fi60W CMD and the Hess diagram for all the stars. The 
most obvious feature in these diagrams is the multi-modal MS that stretches from the turn- 
off (just below saturation) to about three F160W magnitudes below it. Furthermore, the 
MS's breadth and color distribution dramatically changes when moving from bright to faint 
magnitudes. Down to mpieow ~ 20.5, the color distribution across the MSs is consistent 
with what has been previously observed in the optical. As expected from theory, below this, 
the MS runs almost vertically (actually slightly shifting towards bluer colors), and the three 
MS components discovered by P07 appear to merge at this bending point. However, below 
m Fi60W ~ 21, the MS split re-appears, but this time with a different morphology. Whereas 
in the upper part of the CMD, and as in the optical bands, the less-populated MS is the 
bluest component, below mpi6ow ~ 21, the less-populated component is the reddest one. In 
the following, we will refer to these two MS-components fainter than mpi6ow ~ 21, as MSj 
(more-populated and bluer), and MSfi (less-populated MS and redder). 

To estimate the fraction of stars in each MS we followed the procedure illustrated in 
Fig. [2l which is similar to that used in several previous papers (e.g. P07). The left panel 



2 http:www.stsci.eduhstwfc3phot_zp_lbn 
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Fig. 1. — Upper panels: mpnow versus mpnow — mpi 60 w CMD for stars in the fields 'A', 'B', 
and 'C corrected for differential reddening. Lower panels: mpi 60 w versus mpnow — m F160W 
CMD (left) and Hess diagram (right) for all the stars. 

shows a zoom-in of the mp 160 w versus mpnow — mp 160 w CMD region where the MS split 
is most evident. We analyze the color distribution of the stars plotted in the left panel in 
four magnitude intervals over the range 21.25 < mpi6ow < 22.5. The distribution is clearly 
bimodal, and has been fitted with two Gaussians (dashed and continuous black lines for the 
MS\ and MSn, respectively). From the areas under the Gaussians, we find that 65±2% of 
stars belong to the MS\, and 35±2% to the MSn. The errors are calculated as the rms of 
the measurements obtained in the N=4 magnitude intervals divided by \/N — 1. 

We note that the fractions of stars along MS\ and MSn are very similar to the fraction 
of rMS (62±2%) stars and the total fraction of mMS and bMS stars (24+14=38±3%, M12). 
This fact makes it very tempting to associate the MS\ with the rMS (defined by P07), while 
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Fig. 2. — Left panel: A close-up of the region of the mFi60W versus mFiiow — ^Fi60W CMD 
from the lower-left panel of Fig. 1, where the MS split is more evident. Right panel: The 
color distribution of the stars plotted in the left panel, in four F160W magnitude intervals. 
The continuous gray lines are fits by a sum of two Gaussians. 

both their mMS and their bMS could be the extension to bright magnitudes of the MSu. 
To further investigate this issue, in Fig. [3j we show photometry for the stars in field 'A' that 
also happen to be measured in P07. The three groups bMS, mMS, and rMS defined in P07 
are colored in blue, green, and red, and are plotted with the same color-code in the other 
panels. In the upper-right panel, we show the m F160W vs. m F now — m Fi60W CMD for the 
same stars cross-identified in the WFC3/NIR bandpasses. We note that: i) above the MS 
bend (mFi6ow — 20.5) the three MSs exhibit in the NIR the same relative locations as in 
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the optical CMD, ii) At the magnitude of the MS bend, the three MSs appear to merge, 
and hi) Below the MS bend, the bMS and mMS appear to become redder than the stars 
labeled as rMS in the optical photometry. The histograms of the NIR-color-distribution in 
the inset of the top-right panel of Fig. 3 appear to confirm this trend. The cyan and the 
red histograms show the color distribution for bMS+mMS stars and rMS respectively in 
the magnitude interval 20.75< mFi6ow <21.00 The cyan histogram (bMS+mMS) has, on 
average, bluer colors than red histogram (rMS). 

In the lower-left panel we show the NIR fiducial lines for the optically selected rMS 
sample of stars (in red). Defining this fiducial line is an iterative procedure. We started 
from our optically defined sample of rMS stars (from P07 and M12), then we calculated the 
sigma-clipped median of NIR colors and of NIR magnitudes for rMS stars in intervals of 0.2 
mag in F160W, and interpolated with a spline. Similarly, we obtained the fiducial lines for 
the combined sample of mMS and bMS stars shown in the middle panel. The two fiducial 
lines are compared in the lower-right panel and confirm that the rMS is indeed redder than 
the mMS and bMS for magnitudes fainter than mFi6ow ~ 20.6. It is therefore legitimate to 
associate the rMS with the component MSi, and bMS+mMS with the MSn- 



4. Discussion 

The near-IR/WFC3 CMD presented in this Letter allows us to examine, for the first 
time, the behavior of multiple MSs among VLM stars. As expected, in the near-IR, the MS 
stars less massive than ~0.4 M define a sequence with nearly-constant color (see e.g. Baraffe 
et al. 1997, Zoccali et al. 2000, Bono et al. 2010 and references therein). Theoretical models 
predict that this is due to two competing effects. On one hand, the increase of the radiative 
opacity, coupled with the decrease of the effective temperature shifts the stellar colors to the 
red. On the other hand, the increase of the collisional induced absorption (CIA) of the H 2 
molecule in the infrared moves back the stellar flux to the blue. In a given range of stellar 
masses, the two effects compensate, and the MS locus runs almost vertical. Moving towards 
less massive stars the second effect becomes dominant, and the color of the MS becomes 
bluer and bluer with decreasing stellar mass. The CIA source is obviously related to the 
abundance of H2 molecules. It also scales as the square of the density (p), as opposed to 
other opacity sources, which are proportional to p (Cassisi 2011 and references therein). 

In order to compare empirical evidence with suitable evolutionary predictions, we have 
computed some grids of evolutionary models for both low- and very low-mass stars. For 
the low-mass structures (i. e. Ai > 0.5.M©) we adopt the physical scenario described in 
Pietrinferni et al. (2006), while in the VLM regime we use the same physical inputs adopted 
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Fig. 3. — Upper panels: mF8i4w versus mF475w — ^F8i4W (left) and mFi60W versus mpnow — 
^Fi60W CMD for stars in field 'A' (right). The rMS, mMS, and bMS stars are colored red, 
green, and blue, respectively. The color distribution for MS stars between the two dashed- 
dotted lines is shown in the inset where red and cyan histograms correspond to rMS stars 
and to the sample of both mMS and bMS stars respectively (see text for details). Lower- 
panel: Fiducial line for the rMS stars (left) and for stars in the other two MSs (middle) 
superimposed on the IR CMD. A comparison of the two fiducial lines is shown in the right 
panel. 



in Cassisi et al. (2000). We address the interested reader to the quoted references for details. 
The match between the more massive models and the VLM ones was made at a mass 
level where the transition in luminosity and effective temperature between the two regimes 
is smooth (usually ~ O-SAf©). We computed stellar models for an iron content equal to 
[Fe/H] = —1.3, and a— elements enhancement equal to [a/Fe] = +0.4, and for three values 
for the initial He abundances: Y=0.248, 0.35 and 0.40, namely. These initial He abundances 



-9 - 



have been selected on the basis of the comparison between stellar models and the triple MS 
observed in the optical CMD of NGC 2808 performed by P07. 

The theoretical models have been transformed into the observational domain by inte- 
grating the synthetic spectra of the BT-Settl AGSS model atmosphere gricjf] (Allard et al. 
2011, 2012) over the IR WFC3 bandpasses. We used the filter transmission tables provided 
by the System Throughputs web site for the WFC3. 

The upper-left panel of Fig. H] shows the location in the selected near-IR CMD of 12- 
Gyr isochrones for the selected assumptions about the initial He abundances. Data in Fig. H] 
show that for mpi6ow ^ 21 (i. e. M F160W < 6, the He abundance plays a fundamental role in 
driving the MS location: when increasing the He abundance from Y=0.248 to 0.40 the MS 
locus becomes fainter by more than 0.5 mag at fixed color. This is consistent (and indeed 
explains) the separation of the three MSs observed in the upper part (mFi6ow < 20.5) of the 
CMDs of Figs. 1 and 3. However, since helium-rich stars have lower hydrogen abundance, we 
would expect that, in the stellar mass regime where the two - previously quoted - physical 
processes are in competition, the He-rich VLM stars would appear redder than He-normal 
stars for their lower H 2 abundance, and consequent lower contribution coming from the H 2 
CIA. Indeed, at fainter magnitudes, the He-normal (Y=0.248) MS sequence runs on the blue 
side of the more He-rich sequences, but the effect is very small, smaller than the observed 
one. 

We note that, below the faint limit of our photometry (M F160 w J$ 10), the He-normal 
sequence reaches fainter magnitudes and bluer colors than the He-rich counterpart. This 
is an important feature that could be tested with deeper CMDs, and may represent an 
independent observational confirmation of high- He enhancement. 

The bottom-left panel of Fig. H] shows that the difference in the initial He abundances 
among the stars belonging to the distinct MS loci is not able to provide a complete expla- 
nation of the observed trend. In the following, we will attempt to identify the origin of the 
(small) separation of the observed MSs at magnitudes fainter than Mpi6ow ~ 20.5. First of 
all, we note that the model atmospheres adopted for computing the color- T e g- transforma- 
tions do not account for the peculiar chemical patterns of the various sub-populations present 
in NGC 2808. In particular, they do not take into account the effect of light-elements anti- 
correlations, which have been shown to be relevant for the MS color in the UV and Stromgren 
bands (Sll, M12). 



3 These model atmospheres and colors tables are available via the Phoenix web simulator at the following 
URL site: |http: / /phoenix. ens-lyon.fr/simulator 
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Fig. 4. — Left panels: theoretical isochrones for an age of 12 Gyr, a metallicity suitable for 
NGC 2808 and various assumptions about the initial He content (see labels). The upper-left 
panel shows the theoretical absolute magnitudes and colors. In the lower-left panel, we show 
-after applying a distance modulus (m — M) =15.0 and a reddening E(B— V)=0.18, in order 
to make the comparison with the upper right panel of Fig. 3 easier- the location in the CMD 
of the stellar models with mass in the range from O.O894M to 0.65M Q . The arrow marks the 
approximate location along the MS loci of the VLM stellar model with effective temperature 
and gravity values consistent with those adopted for computing the model atmospheres. 
Right panels: Comparison of the synthetic spectra of an bMS star (cyan) and an rMS star 
(red, see text for more details) difference between the synthetic spectrum of a bMS and a 
rMS star (middle). Normalized responses of the F110W and F160W WFC3/NIR HST filters 
(bottom). 

In order to explore the impact of light-element variations on the NIR WFC3 bands, we 
have computed synthetic spectra trying to account for the chemical patterns of the stars 
belonging to the distinct MSs in NGC 2808. For this calculation, we used the ATLAS9 and 
SYNTHE Kurucz programs in the range from 8000A to 18000A (Kurucz 2005, Sbordone et 
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al. 2007)0. 

For the MS\ and MSu, we adopted a Helium content of Y=0.25 and Y=0.38 respec- 
tively, as suggested by isochrone-fitting on the upper MS (M12). For the MS\, we assumed 
average chemical abundance of O-rich stars ([O/Fe]=0.4, as measured by Carretta et al. 
2006), and adopted the carbon and nitrogen abundance ([C/Fe]=— 0.3, [N/Fe]=0.5) mea- 
sured by Bragaglia et al. (2010). For the MSu, we used [0/Fe]=— 0.6 (the O abundance 
measured by Carretta et al. 2006 for the O-poor group), and [C/Fe]=— 0.7, and [N/Fe]=2.0 
(as measured by Bragaglia et al. 2010 for a bMS star). Table 1 summarizes the adopted chem- 
ical abundances for the two sequences. For both MSs we used T c g = 4000K, log(g)=5.0, and 
a microturbolence 1.0 km s _1 . Our synthesis includes the following molecules in the Kurucz 
compilation: CO, C 2 , CN, OH, MgH, SiH, H 2 0, TiO (H 2 0, from Partridge & Schwenke 
1997; TiO from Schwenke 1998) VO and ZrO (B. Plez priv. communication) The result- 
ing synthetic spectra (Fig. 4, upper-right panel) have been integrated over the transmission 
of the WFC3/IR F110W and the F160W filters (lower-right panel) to produce synthetic 
magnitudes and colors. We found a A(Mfhow — ^Fi60w) = 0.10 color difference between 
the simulated MSj and MSu stars, consistent with the observed color difference (Fig. 2) at 
mFi6ow=21-5: A(m F110 w — ?Ti Fi60w) = 0.06±0.01 mag. The middle-right panel of Fig. 4 shows 
the flux difference as a function of the wavelength between the two MSs spectra. The H 2 
molecules have the strongest effect on the synthetic spectra, and cause the significant lower 
flux of the rMS in the filter F160W. We consider the agreement between the simulated and 
observed color differences satisfactory, accounting for the high sensitivity to light-elements 
abundances of the stellar spectra at these wavelengths, and that the adopted C and N abun- 
dances are based on the measurement on just one bMS and one rMS star (Bragaglia et al. 
2010). 



4 http:/ /wwwuser. oat. ts.astro.it/castelli/ 



MS 


Y 


[C/Fe] 


[N/Fe] 


[O/Fe] 


rMS 


0.25 


-0.3 


0.5 


0.4 


rnMS 


0.32 


NA 


NA 


0.0 


bMS 


0.38 


-0.7 


2.0 


-0.6 


MSI 


0.25 


-0.3 


0.5 


0.4 


MSII 


0.38 


-0.7 


2.0 


-0.6 



Table 1: Average chemical abundances of bMS, mMS, and rMS stars and abundanced 
adopted for the MS\ and the MSu- 
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5. Summary 

The photometry presented in this paper for NGC 2808 reveals, for the first time, multiple 
sequences in near-infrared CMD of a GC. The brightest part of the CMD (mFi6ow < 20.5) is 
consistent with three populations with different He and light-element abundance, as already 
noticed in previous papers based on visual photometry. The three MSs merge together at the 
luminosity of the MS bend while at fainter magnitudes, a combination of stellar structure 
and atmospheric effects makes the distinction among the different stellar populations more 
intricate. Our CMD allows us to identify at least two MSs. A redder, more populated 
MSj, which includes ~ 65% of the MS stars, and that we associate with the first stellar 
generation, which has primordial He, and O-C-rich/N-poor stars, and a MSn, with ~ 35% 
of stars, corresponding to a second generation stellar population that is enriched in He and 
N and depleted in C and O. The MS\ of Fig. 1 is the faint counterpart of the rMS identified 
by P07, whereas the MSn corresponds to the lower mass counterpart of the mMS and bMS 
of P07. 

This Letter provides the first detection of multiple populations with different helium 
and light-element abundances among VLM stars and extends the investigation by SI 1 to the 
near-IR. The fact that the signatures of abundance anticorrelation are also observed among 
fully-convective M-dwarfs demonstrates, once and for all, that they have primordial origin 
and hence correspond to different stellar generations. 
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